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ABSTRACT 



The pressure and phase distribution of sound in a fast fluid 
medium underlying a tapered fluid medium was modeled using a Green's 
function approach. The model predicted well defined beams in the 
bottom, with apparent interference effects strongly evident for cases 
at a distance from the wedge interface. The effects of attenuation on 
the patterns were studied. A simple expression for determining the 
beam depression angle as a function the wedge angle, the sound speed 
ratio, and the density ratio was derived. Comparisons with laboratory 
measurements were made. 
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I . INTRODUCTION 



The problem of radiation behavior in a wedge-shaped medium 
overlying a bottom possesing a greater speed of sound has been 
investigated both optically and acoustically. In 1971, Tien and 
Martin (Ref. 1) deposited a thin, tapered, dielectric film on a 
substrate having a higher refractive index, and observed the behavior 
of a laser beam coupled into the film. Perfect reflection at the 
film- substrate interface was observed until the changing angle of 
incidence decreased below the critical angle for continued mode 
propagation. For angles of incidence below critical, the light was 
converted into radiation modes in the substrate. 

A concurrent analysis of the acoustic problem was done by Kuznetzov 
(Ref. 2). Kuznetzov developed a theory, based on normal-modes, that 
predicted the following for sound propagating within a wedge-shaped 
medium overlying a half-space with a higher soxind speed: 

1. As the sound approached the vertex of the wedge, it would 
continue to be totally reflected until the angle of incidence decreased 
to the limiting angle of total reflection. 

2. Sound incident at less than the limiting angle of total 
reflection would be totally refracted into the underlying half-space. 
This refraction occurs within the region from the wedge apex to the 
point where the limiting angle was attained. 
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3. Within the half-space , the acoustic energy would he propagated 
as a well defined beam with the maximum pressure amplitude occurring at 
an angle of depression from the wedge interface within the range 3 to 23/ 
where 3 is the wedge angle. 

Kuznetzov also performed a series of experiments that supported his 
theory . 

Subsequent optical experiments by Tien, Smolinsky , and Martin (Ref. 3) 
demonstrated the formation of a well-defined beam within the substrate . 

The accompanying theoretical discussion predicted conflicting behavior . 
Ray-optics predicted refraction into the substrate after the cutoff dis- 
tance (the point where refraction begins) and wave theory predicted com- 
plete refraction into the substrate before the cutoff point. Sigelman , 

et. al. (Ref. 4) both predicted and observed a well columnated beam in a 

water- aluminum wedge system, but the observed beam was much broader than 
predicted. Also, the pressure maxima for the columnated beam occurred at 
approximately 11^ for a wedge angle of 1.3^^, well outside the limit pre- 
dicted be Kuznetzov. 

The primary objective of this theoretical development is to produce 
a simple, fast computer model to predict the pressure and phase distri- 
bution of the columnated sound beam in the fast substrate. The pressure 
and phase distribution along the wedge interface calculated by Kawamura 
(Ref. 5) will be used as input to the substrate model, and will specifi- 
cally be used as the boundary conditions for a Green's function analysis . 

To remove explicit frequency dependence, all relevant distances will be 
normalized to the distance measured from the apex along the wedge inter- 
face at which the critical reflection angle is first exceeded; this dump 
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distance is referred to as the dump distance X. Integration over the 
wedge interface will be truncated at distances sufficient to encompass 
only selected modes of interest . Results will be analyzed for trends 
and compared with existing measured data. 
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II. THEORY 



Pressure and phase distributions along the bottom of a wedge- 
shaped medium are obtained using an adaptation of the method of images 
as prescribed in Ref. 5. The source may be set at either finite or 
infinite distance, but must be maintained sufficiently far from the 
wedge boundaries to ensure that impinging waves are essentially locally 
planar. The following theoretical discussion uses the above pressure 
and phase distribution as initial conditions, and calculates the 
pressure and phase distribution in the fast bottom by utilizing the 
Green’s function (Ref. 6). 

Figure II-l illustrates the geometry of the problem. The horizontal 
x-axis and the vertical z-axis have their point of origin at the wedge 
apex, with the positive x-axis forming the interface between the wedge 
and fast bottom. The wedge angle is designated B and the angle 
specifying the sound source location is 6 , X is the first mode 
dumping distance, and is a distance determined as follows: 



0 



arccos ( C 2 / Ci ) 



(Il-la) 



c 



D 



0 



Cl / { 4*f«sin(0^J } 



(Il-lb) 



X 



/ sin(B) 



(II-lc) 
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Figure II- 1 Wedge Problem Geometry 



where 



Cl = sound speed in the wedge 

C 2 = sound speed in the bottom 

f = frequency 

and = lowest mode cut-off depth. 

To make the problem frequency independent, all distances were 
normalized to X by defining 

C = X / X (II-2a) 

r = xW X (II-2b) 

and n = z / X Cll-2c) 

The variable is identical to x, and is identical to The 
prime designation is required for dummy variables of integration in 
the Green's function development to follow. 

Pressure along the wedge bottom is then given by 

p(5,n=0,t) = P(5) exp {i[ojt + L(C)1} (II-3) 



where P(5) is the pressure amplitude along the wedge bottom and L(55 
is the phase along the wedge bottom. Pressure in the bottom, away 
from the boundary, is then given by Ref, 7, 



/ 

p(r,t) = exp[io3t] ^ j P(0 



8G 

dz^ 



X d5^ 



z"=0 



= exp [io)t] R(r) 



(II-4) 



11 



where G is the appropriate Green’s function. The z* coordinate is 
identical to the z coordinate and r = (x^ + The appropriate 

Green's function, for a line source in cylindrical geometry is 

G(r,r") = (k 2 1 r - r" | ) (II-5) 

where J <2 is the propagation vector in the bottom and r = [ (x' ) ^ + (z ' ) . 

Differentiation of the Green's function yields (Ref. 7) 

, V .V k 2 (z - z") 

0^2 ht “ r" I ) (II-6) 

1 -^1 ( i ) 

where k 2 | r - r | >> 2ir and , H, are zero and first order Bessel 
functions of the third kind (Hankel functions) respectively. Using 
equations (II-2a,b,c) at z^ = 0 the expression [r - | may be 

rewritten as 



r - r 



= X [ (C - 5^)2 + 



Cll-7) 



3 G 

Substituting (II-7) into equation (II-6) ^ and evaluating at z'' = 0 

d Z ' 



yields 



9G 

3z" 



'z"=0 



^ + n^]^}k2nX 



(II-8) 



Combining equations (II-8) and (11-4) yields a time independent 



expression for amplitude and phase, R(r) 

. . . 1/2 

R(r) = ^ exp 






P (Z' ) exp{j[L(^' )-k:>X 



(II -9) 



Equation (II-9) is the basic equation used for the calculation of 
phase and amplitude in the bottom. 
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Ill, DEVELOPMENT OF WDGBTM 



Program WDGBTM was designed to perform four main functions . 

First, the pressure and phase along the interface of a wedge and 
underlying fast media may be calculated in two ways, for both 
infinitely and finitely distant sound sources in accordance with 
Ref. 5. Modifications were made to the code of Ref. 5 to facilitate 
subsequent in--bottom calculations, and to enhance user flexibility. 

These changes did not affect results . Then using the first code as 
input, additional code enables both the computation of pressure and 
phase along a line perpendicular to the wedge bottom, and as a function 
of radius and angle measured from the first mode dumping point, X. 

Modes of program execution may be chosen by literal input directives 
combined with appropriate numerical input. Appendix A. illustrates an 
input set that utilizes all four options. Output consists of both 
printed listings and graphics. 

Program versions exist for execution on both a CDC 6500 and an IBM 
360. The CDC 6500, with a 60 bit word, does not require double 
precision for accuracy equivalent to the IBM 360 output. The CDC 
version of WDGBTM has simpler varian control routines, and separate 
output files for each pressure/phase variation# This enables optional 
suppression of listings. Execution time on the CDC 6500 is roughly ten 
times faster than on the IBM 360. For average input sets, times are 150 
and 1500 cp seconds respectively. 
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A. EQUATION PREPARATION 



1. Depth Domain 

Integrations along lines perpendicular to the wedge bottom 
were designated depth domain calculations. The program was designed to 
allow specification of the field point, i.e., specification of where 
along the media interface line the perpendicular probe line could be 
dropped. The field point is not confined to lines directly beneath the 
wedge, but may extend to points beyond the wedge apex. These points 
would be values of negative ^ in Figure II-l. The beginning and ending 
points of the section of the perpendicular to be examined, and the 
mamber of field points used in the interval were also to be specified. 

Equation (II-9) was split into real and imaginary components before 
coding. Integration order selected was over the entire C ^domain for 
each n point. 

2 . Angle Domain 

Examination of Figure ir-l verifies the following geometric 
relationships. 



and 



€ = O cosQ, 
k 



(Ill-la) 



n = iiii-lb) 

where £ is the distance from X to the coordinate of the field point 
of interest and a is the radial distance from X to the field point. 
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Substitution of equations (Ill-la) and (Ill-lb) into equation (II-9) 
yields 



i — [- 4 -) 



(III-2) 



r? (V )exp{i[L(V ) ~k2X[(E,~V 4 +0^ -Z^V^^]}dV 

Equation (III-2) is the basis for angle domain calculations. The 
specification of a originating at X was chosen to conform to existing 
experimental data (Reference 8) . The length of the radius, the 
beginning and ending of the angular sweep, and the angular increment 
comprise the input specifications. As in the depth domain option, 
equation (III-2) was split into real and imaginary components before 
coding. 



3 . Design Restrictions 

Theoretical development of program WDGBTM assumes the sound 
source is sufficiently distant to ensure essentially locally planar 
waves at the interface and the in-bottom measurements satisfy the 
condition Ir-r'^>>2TT, 

4. Program Results 

An initial set of WDGBTM runs were made in the depth domain for 
a perpendicular directly down from the wedge apex. The source was 
infinitely distant and set at an angle of 6/2. This limited modes within 
the wedge to odd numbered (symmetric) modes only. Integration over the 
wedge bottom was truncated at C equal four, so that only modes one and 
three were considered as contributing significantly to the beam pattern 
in the bottom. Beam patterns in the bottom were calculated over the range 
n equal zero to 1.3X. 
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The initial set was comprised of seventy- five runs. In addition 
to the above criteria / the runs were generated by all possible combi- 
nations of the following specifications. The ration P 2 0.95, 

0.90, 0.80, 0.70 or 0.50. The sound speed ratio, 0.95, 

0.90, or 0.85. The wedge angle 6 was varied through 10, 6.92, 5, 2 and 
1 degree (s) . 

From the results of the initial data set, the depression angle and 
beam width was calculated for the beam pattern in the bottom. Results 
are summarized in Appendix B. These results were then used to find a 
simple analytical formula to predict the depression angle as discussed 
in Section IV. 

Other runs were made for perpendiculars at x - -9X from the wedge 
apex. These runs displayed apparent strong interference phenomena and 
are discussed in detail later. 

Angle domain executions were made to compare with experimental data. 
Results are discussed in Section VII. 

B. SENSITIVITY TO INTEGRATION STEP SIZE 

It was desirable to keep both the integration step size and the field 
point increment size as large as possible in order to minimize total 
computer resources required. Increment of 0.01 times the dump distance 
(X) proved to be optimum for both integration increments along the 
wedge interface and field point steps in the bottom. 

To illustrate the adequacy of O.OIX increment, a series of runs were 
made in accordance with Table III-l. An overlay of wedge interface 
pressure amplitude for = O.OIX and Ax = O.OOIX demonstrated no dis- 
tinguishable difference. Likewise , depth domain plots (a, b, c, d of 
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Table (III - 1) 



Sensitivity to Integration Step Size (A^) and Field point interval 



size (Ati) 

A;f = 0.001 




An=o.oi An=o.oo2 



Ax = 0.0001 




An=o . 01 An=d . 001 



(a) (b) 



(c) (d) 



Run 


T1 amp max 


CD 

a 


a 


0.240 


13.50 


b 


0.244 


13.71 


c ' 


0.240 


13.50 


d 


0.243 


13.66 
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Table could be identically superimposed. Use of an Integration 

step size of 0,001X resulted In a calculated depression angle within 
one percent of the depression angle for the O.OIX run. Utilization of 
a step size of greater than 0,01X Is not recommended. Within this 
limitation the Integrations can be carried out over lOX along the Inter- 
face for a maximum field Interval of lOX In the bottom. The lOX re- 
striction Is Imposed by program array size, (For the 0^'S specified ^ 0,5 
degrees In the angle domain was slightly finer than the 0,01X field 
point step size In the depth domain and was the field point Increment of 
choice) , 
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IV. BEAM DEPRESSION ANGLE ANALYSIS 



A simple relationship was sought between the physical parameters 
associated with the wedge system and the resultant beam. 



Assume 




• 


0OC . B” 



where 9„ is the beam depression angle and B is the wedge angle. For 

LA^ 

each P^/P 2 ^ graphed on log-log paper, and the 



slope ( n) was measured. 


Results are listed in Table 1 of Appendix C, 


n equals 0.329. 





Then found as a function of The critical angle, 9^, 

is an explicit function of 

9^ = arccos (C^/C^) ,9^ in radians (IV-2) 



It was then assumed 






^DC "" 



For any constant B , M was derived as the slope of a log-log graph of 
®rc Graphs were analysed for 6 equal to 10, 6.92, 5, 2 and 1 

degree, and 0.95, 0.90, 0.85, 0.70 and 0.50. Results 

are summarized in Table (2) of Appendix C, ^ - 0.772. 
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The same process was repeated for ascertaining dependence on 
P 1 /P 2 • Log- log graphs were prepared by holding B and 9^ constant 
using the assumed relation 



'DC 



e. 






(IV-4) 



Results are summarized in Table 3 of Appendix C. q - -0.254. 

The proposed formula is now dependent on all physical properties, 
assuming the source angle is maintained at one-half B. To remove the 
proportionality sign, write 



and solve for K , 
Therefore 



0DC = ^ (iv-5) 



From Table 4, Appendix C, 17.22/radian 



. 17.22 (iv-6) 

where B is in degrees and radians. 

0^ calculated from equation (IV-6) and the difference between 9^^ 

and 0 from program WDGBTM (labeled A0 and 9 respectively) are also 
p D DW 

tabulated on Table 4, Appendix C. A0^ varies from a minimum of 0.00 to 
a maximum of 1.25, with the average of the absolute values of A0^ equal 
to 0.36. Errors ranged from zero percent to 5.0 percent. 2hus ^ gives 
a reasonable approximation of the beam depression angle. 
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A similar formula 



e,, - I..I a/-- rvp,r-- 



{lV-7) 



was developed independently by A. Coppens (Reference 6) The minimum 

deviation of from 9^ is 0.01 and the maximum deviation is 1.56. 

DA DW 

The average of the absolute value of the deviation is 0.36. Both 
formulas are best for P 2 close to 1.0, and less accurate for 
p^/p^ equal to 0.5. 
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V. THE EFFECT OF ATTENUATION 



As previewed in Section III-A.4 the pressure distribution as a 
function of U at ^ = -9.0X demonstrated apparent interference effects 
(see figure V-1) . As integration over the wedge interface to 4X 
encompasses dumping by the first and third mode, and the third mode 
diimps about the point 3X vice X for the first mode, it is logical to 
expect attenuation effects to reduce the contribution of the more 
distant mode. The initial Green’s function approach, equation (II-9) 
did not consider attenuation of the acoustic wave in the bottom 
medium. An attenuation term 

exp (v-1) 

was incorporated into equation (II-9) to yield 

1/2 

R(r) = r\/4 exp^^'j (V-2) 

f 'p(^')exp{-oi’ )^+r\^^^^}exp{i[L(Z' ) -k2X 

An appropriate code modification to the depth domain routine was made. 
The attenuation coefficient ct was an input variable to the program. 

Two sets of comparison riins were made. The first set consisted of 
an unattenuated control (0t= 0) and an attenuated run with a equal 1.0, 
with the pressure amplitude measured down from the wedge apex. The 
second set consisted of an unattenuated control plus runs with a set to 
1.0, 0.5, 0.1, and 0.01 and the pressure measured at = -9.0 from the 
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apex. All other input variables for both sets were held constant and 
are specified in Table (V-1) . Integration over the wedge bottom was 
from the apex to four times the dump distance; this was to insure that 
contributions from the third mode (dumping about the point x- 3X) were 
included. 

The results of set one are illustrated in Figure (V-2) . Note that 
the second maxima is severely depressed in the attenuated case. The 
depression angle has been shifted slightly from 11.86 degrees for the 
unattenuated data to 11.31 degrees with ct set to 1.0. 

Figure V-3 is a plot of the results from set two. All curves have 
been normalized to the initial slope for Ct = 1.0. The second and third 
maxima decrease with increasing Ot, and relatively more energy is 
contained in the first beam. The depression angle of the first beam 
decreases slightly with increasing a . These trends are summarized in 
Table (V-2). 

Figures (V-4a) and (V-4b) illustrate the effects of attenuation on 
different contributing modes. For Figure (V-4a) , the source was placed 
at 1/3 S. With this source geometry, contributions from mode three are 
suppressed, and only the first and second mode contributed to the beam 
pattern in the bottom. For Figure (V-4b) , the source was placed at 
and modes one and three propagate into the bottom medium. Both cases 
were run for a = 1.0, 0.5, 0.1 and 0.0. As could be expected, 
attenuation causes greater suppression of third mode interference 
contributions than of second mode contributions. 
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Table (V-1) 



Input Specifications for Attenuation Runs 



= 1500.0 
= 1578.947 

2 

C^/C^ =0.95 
6 = 5.0” 
X2 = 4.0 

6 = 2.5“’ 



= 1.0 

= 1.052632 
P^/P^ = 0.95 

f = 150000 hz 
Source distance = 



a = 0.0, 0.01, 0.1, 0.5, 1.0 



Refer to Figure II-l for variable meaning. 



infinite 
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Figure V-2 



Attenuated and unattenuated pressure amplitudes 
down from the apex 
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Figure V-3 



OL = 0 
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TABLE V-2 



ot 


Y 




h 


^2 




1.0 


12.84 




0.21xl0~‘^ 


NONE 




0.5 


13.39 




.33xl0~^ 


.26xl0~^ 


1.26 


0.1 


13.82 




0.21 


0.23 


0.90 


0.01 


13.93 




0.54 


0.66 


0.83 


0.0 


13.98 




0.61 


0.74 


0.82 


Y = 


depression angle of 


first maxima 






h = 


amplitude of 


first maxima 






II 


amplitude of 


second 


maxima 
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BRMP 



RHOl 


0.9651 


EBEG 


0.0 


RH02 


1 . 1470 


EENO 


4.0 


RRflT 


0.8415 


EINC 


400. 


Cl 


1500.0000 


FPNT 


o 

o 

o 


C2 


-6^4. 1984 






CRRT 


0.8854 


Figure 


v-4a 


BETA 


1.53 






SRCfl 


n.51 






SRCO 


999.00 






HINC 


O.OiO 


a 








0.0 





0.1 



0.5 



1.0 
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BflMP 



RHOl 


0.9651 


EBEG 


0.0 


RH02 


1 . m ?o 


EENO 


4.0 


RRflT 


0.8415 


ElNC 


400. 


Cl 


1500. OOOQ 


FPNT 


0.00 


C2 


1694. 1984 


Figure 


V-4b 


CRRT 


0.8854 






aCTFi 


1.53 






sacR 


0.76 






SRCD 


999.00 






WINC 


0.010 


a 





0.0 



0.1 



0.5 



1.0 
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VI. ANALYSES OF FAR-FIELD INTERFERENCE PATTERNS 



The results of the attenuation investigation were supportive of an 
interference explanation for the multiple beams. Subsequently, two 
analytical approaches were tried to predict the location of interference 
maxima and minima under unattenuated far-field conditions. 

Both the Differential Change and Multiple Lambda analyses assume 
very simple processes. Both assume point sources of sound at the wedge 
interface located at X and 3X. Reallistically the sound sources behave 
as phased arrays propagating assymetrically about the dumping points. 

Phase differences between the two dumping points are non-zero. The 
Differential Change formiula becomes invalid as f) approaches zero. 

A. DIFFERENTIAL CHANGE APPROACH 

Based on the simple geometric considerations illustrated in Fig (VI-1), 
the following nondimensional formula for the change in interference 
maxima locations was derived. 



all other variables held constant. Table (VI— 1) is a summary of one case, 
this case is typical of results observed. Results will be discussed in 
Section VI-C. 






Note that equation (VI-1) is dependent only on tne wedge angle p and the 




sound speeds of the two fluid media. 



A series of WDGBTM executions were made for varying and S, with 
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Source at 1/2 
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TABLE VI-1 



Comparison of An from equation (VI-1) and An from WDGBTM results. 

^i/Ci = 0.90 



6 


WDGBTM 


An 


An 




maxima 

locations 


from 

WDGBTM 


predict 


5.0 


2.38 


1.67 


4.56 




4.05 

5.97 


1.92 


3.04 


2.5 


2.33 


1.40 


2.33 




3.73 


1.02 


1.60 




4.75 
5 . 66 


0. 91 


1.38 


1.25 


2.08 


1.22 


1.29 




3.30 


0.66 


0.88 




3.96 

4.45 


0.49 


0.77 
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B. MULTIPLE LAMBDA APPROACH 



The second interference pattern calculation also utilizes the 
geometry of Fig, (VI-1). For each beta combination, 

X and (the wavelength in the bottom media) was calculated according 
to the equations 



4f sin 9^ sin S 

and 

Xj - c/f 

where is ISOkhz. The ratio X^/X = Q 
is then found and used in the relation 

NQ = (vi-4) 

where ^ is any positive integer. Note that use of equation (VI-4) 
assumes zero initial phase. 

To find the expected maxima and minima, a plot of H vs ^ was made 
for each Q. Maxima are expected for each integer iV, minima for each 
^ +1/2. Table (VI-2) summarizes one case of maxima and minima expected 
as a result of equation (VI-4) versus maxima and minima observed on 
WDGBTM output. 

C. SUMMARY OF ANALYSES 

Predictions from both simple analyses are of the same order of 
magnitude as WDGBTM results. The number of maxima and minima are 



(VI-2) 



(VI-3) 
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TABLE VI- 2 



OCCURRENCE OF MAXIMA AND MINIMA 







C 1 /C 2 =-0.85 






^(DEG) 


n MAX WDGBTM 


n MAX EXP 


n MIN WDGBTM 


n MIN EXP 


5 


4.03 


2.68 


5.39 


4.70 






6.34 






2.5 


3.59 


2.59 


4.55 


3.76 




5.00 


4.70 


5.66 


5.56 


1.25 


2.02 


2.53 


2.49 


1.76 




2.85 




3.56 


3.16 




3.81 


3.72 


4.77 


4.25 




5.03 


4.68 


5.40 


5.13 




5.62 


5.54 


5.90 


5.94 
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consistantly predicted as increasing with decreasing wedge angle. 
Occasinal unexplainable phenomena, such as increasing WDGBTM At) for 
the S = 5 on Table (VI-1) and the extra WDGBTM maxima for 3=5 
Table (VI-2) , were observed, but results are generally supportive of 
an interference effect occurring. 
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VII. COMPARISONS 



A. EXPERIMENTAL DATA 

Netzorg (Reference 4) measured pressure amplitude as a function of 
angle for fixed radii in four media below a wedge shaped upper layer. 

The angle domain version of WDGBTM was used to compare calculated 
pressure patterns with measured patterns. 

All runs utilized finite source input. For each of Netzorg *s four 
experimental cases, corresponding computer runs were made with the 
source distance set to lOOX, lOX and SX, where S corresponds to the 
actual source distance used in the tank. In all cases S was less than 
lOX. The integrations were carried out over one dumping distance and 
four dumping distances for each source. 

Figure (VII-1) is typical of pressure amplitudes in the bottom. 

The source in the experiment was only about 40 wavelengths from the 
dumping point, even at lOX the source is approximately 110 wavelengths 
away. Since the method of images calculations utilized to determine 
the pressure amplitude along the wedge bottom assume a plane wave, it 
is questionable that this program can produce a realistic comparison 
with the measured data. 

Figures (VII-2) and (VII-3) illustrate the wide variation in results 
obtained depending on the source distance and integration interval. 
Calculations demonstrate expected trends, i.e., pressure amplitude 
increases with decreasing radius distance (C) , with closing source 
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distance, and as a function of an increasing integration interval over 
the wedge bottom. The beam depression angle becomes more shallow as 
source distance increases. Results at very low angles are not 
considered valid because of assumptions inherent in the WDGBTM model. 

The circled data points of Figures (VII-2) and (VII-3) are a 
transposition of data from Reference 4, Since absolute amplitude was 
not specified in the reference, rescaling to the case most near to the 
actual conditions was performed, i.e., the data were normalized to the 
case = 1,0 and srcd = 3.5X. The shape of the experimental curve and 
the apparent beam depression angle are the only valid comparisons that 
should be made. 

Comparisons of WDGBTM beam depression angles with measured values 
gave varying results. Measured depression angles decreased with 
increasing O , WDGBTM depression angles increased. WDGBTM results were 
within 6 percent to 62 percent of measured values. No correlation was 
observed between prediction accuracy and source distance or prediction 
accuracy and the interface integration interval. All measured beam 
depression angles in Reference S were for the two shortest radii of 
each case, the three longer radii could only be qualitatively compared 
with graphical data. Possible attenuation effects were not considered. 
Numerical results are contained in Appendix D, 

B, PARABOLIC EQUATION APPROACH 

Kuperman (Reference 9) describes a Parabolic Equation model result 
for a wedge problem, and compares these results to one set of data 
extracted from Reference 10. 
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The Parabolic Equation results show very well defined beams in the 
bottom, with no apparent interference effects as observed in comparable 
WDGBTM riins (Figures V-1 and V-4) . Kuperman used an attenuation factor 
of O.SdS per wave length, a value equivalent to an input of 0.984 
Neper s/X for WDGBTM. It should be noted in Figure V-4 that there is 
complete separation of the beams for WDGBTM when ct = 1.0. 

Kuperman states agreement within 20 percent of a measured data set 
of Reference 10. This measured data is identical to Case 1 of Appendix 
D. WDGBTM accuracy for this case varied from zero percent error to 22 
percent error with error calculated according to the formula 

(obseirved data) - (measured data) 
measured data 
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VIII. CONCLUSIONS 



1. Program WDGBTM offers a fast, efficient means of analysing the 
behavior of sound in a fast bottom underlying a wedge-shaped medium. 

Effects of attenuation are easily demonstrated, and modes desired are 
readily controlled by specification of source location and integration 
interval selected. Results of the model are not inconsistant with 
current measured data, and illuminate observations generated by the 
Parabolic Equation model. 

2. A simple equation for calculating beam depression angle as a 
function of wedge characteristics was obtained which offers sufficient 
accuracy for use as a rapid analytical tool. The equation was derived 

for acoustic sources placed at one-half the wedge angle. Accuracy degrades 
when the wedge medium density to bottom medium density ratio decreases 
to 0.5. 

3. It is recommended that further experiments be performed that 
would further test the model. New experiments should be sufficiently 

large scale to ensure plane waves at the wedge bottom interface and to satis- 
fy the k^\r-r |>>27r criteria. Bottom media that reasonably model 
real ocean fluid-like bottom are also recommended. 
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APPENDIX A 



SAMPLE INPUT SET 





O 


o 






o 


o 


o 


o 

• 

\T\ 


>c 




-o 

• 

OvJ 



o 

m 



iTk 




















ro 










>o 










CNJ 










in 






00 fO 




om 






om 




• • 






• • 
















O 


o 








• 


• 




o 




(T' 


O' 




• 




1 


1 




o 



o • 


o • 






oo 


oo 






oo 


oo 






oo 


QOO 






oo 


0^0 


in 


in 


• ITS 


•in 


• 


« 


1-^ 1—4 


O f— < 


o 


o 



r^ 








o 


o 


O 


0 

0 


1^ 


o 


O 


• • o 




• 


• 


in in • 




o 


o o 


>i- v*H 


— \ • 




• 




CLOD 




H-QO 




zr^o 


»— 


H- OCVJO 


H- K H- 


►-•in • 


o 


Z) Q.O • 


o o z> 




CL 


CL 2--<^ 


a <x CL 


tu 


2 


2 •-• 


2 2 2 


u 


»— « 


►>H 


NH NN 


cc 




tu 


O O 


o 


2 


2 O 


2 *2 #2 


o 




•-• QC 


hhO*—» O*-^ 


l/> 


< 


< z> 


< < < 




2 


2 O 


2 2 2 


tu 


a 


o 00 


o o o 


»- 


Q 


Q 


O Q Q 


— lO 




LUO 


O 


2 • 


X 


OX o^- • 


• tu tu X 




OH- 


Oh~ c 


j < 



uL000CL0»^Q-00^000r^0m0c>a.00 
ZlTk *>t*LU • •tu • •«— »m •LU • • 

QOOQOOu.'^O <*^<0000 



c; 



*s 

o 






c; 

o 



44 



Litoral string for infinite source 
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APPENDIX B 



P^/P^ = .9500 



c^/c^ 



6 (DEG) 



DEPRESSION 



BANDWIDTH 



RATIO 



0.95 
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= ARCTAN(Z/X) AT MAX 
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AKPLITUDE 





BANDWIDTH 

- 4-d 

MEAN BW 



= AECTAN(Z/Xj^^ _ 
= .73 
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9000 



c^/c^ 
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APPENDIX C 



TABLE 1 



^ 1/^2 


C /C 

r 2 




n 


0.95 


0.95 




0.318 


0.95 


0.90 




0.347 


0.95 


0.85 




0.341 


0.90 


0.95 




0.318 


0.90 


0.90 




0.335 


0.90 


0.85 




0.339 


0.80 


0.95 




0.328 


0.80 


0.90 




0.328 
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0.50 


0.95 




0.339 


0.50 


0.90 




0.314 


0.50 


0.85 




0.314 






Average n 


= 0.329 




Standard (N-1) 


deviation 


= 0.014 



(N-1) Variance = 0.002 



51 



TABLE 2 



^1/^2 


6 


M 


0.95 


10.0 


0.834 


0.95 


6.92 


0.845 


0.95 


5.00 


0.783 


0.95 


2.00 


0.739 


0.95 


1.00 


0.752 


0.90 


10.0 


0.815 


0.90 


6.92 
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0.90 
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0.90 
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0.90 
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0.80 
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0.80 


5.00 
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0.80 


2.00 
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0.80 
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0.750 


0.70 


10.0 
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0.70 


6.92 


0.892 


0.70 


5.00 


0.701 


0.70 


2.00 


0.778 


0.70 


1.00 


0.697 


0.50 


10.0 


0.786 


0.50 


6.92 


0.881 


0.50 


5.00 


0.693 


0.50 


2.00 


1.051* 


0.50 


1.00 


0.722 


Rejected data point Average M 


= 0.772 




Standard (N-1) deviation 


= 0.059 




Variance 


= 0.003 
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Table 3 





0^(rad) 




B 




0.95 


0.318 




10.0 


-0.217 


0.95 


0.318 




6.92 


-0.285 


0.95 


0.318 




5.00 


-0.257 


0.95 


0.318 




2.00 


-0.215 


0.95 


0.318 




1.00 


-0.283 


0.90 


0.451 




10.0 


-0.183 


0.90 


0.451 




6.92 


-0.232 


0.90 


0.451 




5.00 


-0.211 


0.90 


0.451 




2.00 


-0.300 


0.90 


0.451 




1.00 


-0.317 


0.85 


0.555 




10.0 


-0. 162 


0.85 


0.555 




6.92 


-0.242 


0.85 


0.555 




5.00 


-0.255 


0.85 


0.555 




2.00 


-0.362 


0.85 


0.555 




1.00 


-0.283 






Average 7 


= -0.254 






standard (N-1) 


deviation 


= 0.053 








variation 


= 0.0028 
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Table 4 



1/K 2 






®BW 


K 


®DC 


A0d 


0.95 


0.95 


10.0 


15.64 


17.54 


15.35 


- 0 .29 


0.95 


0 . 95 


6.92 


13.50 


17.08 


13.60 


0.10 


0.95 


0.95 


5.00 


12 .02 


16 . 94 
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0.20 


0.95 


0.95 


2.00 


9.09 


17.32 


9.04 


- 0.05 


0.95 


0.95 


1.00 


7.58 


18.14 


7.20 


- 0.38 


0.95 


0.90 


10.0 


19 .80 


16 .94 
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0.95 


0.90 


6.92 


17.74 


17.13 


17.83 


0.09 


0.95 


0.90 
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16 .33 


17.55 
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0.95 


0.90 


2 . 00 
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0.95 


0.90 
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16.59 


9.43 


0.32 
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6 . 92 


21.16 


17 . 42 


20.92 


- 0,24 


0.95 


0 .85 


5.00 


18.42 


16.87 


18.80 


0.38 


0.95 


0 . 85 


2.00 


13 . 50 


16.72 


13 . 91 


0.41 


0.95 


0 .85 


1.00 


11.31 


17.59 


11.07 


- 0,24 


0.90 


0.95 


10.0 


15.64 


17.31 


15.56 


- 0.08 


0 . 90 


0.95 


6.92 


13.49 


16 .85 


13.79 


0.30 


0.90 


0.95 


5 .00 


12.02 


16.71 


12.39 


0.37 


0 .90 


0.95 


2 .00 


9.09 


17.08 


9.16 


0.07 


0 . 90 


0 .95 


1.00 


7.58 


17.89 


7.30 


- 0.28 


0.90 


0.90 


10 . 0 


19.80 


16.71 


20.40 


0.60 


0.90 


0.90 


6.92 


18.42 


17.55 


18.08 


- 0.34 


0 .90 


0.90 


5.00 


16.33 


17.31 


16.24 


- 0.09 


0.90 


0.90 


2.00 


12.02 


17.23 


12.02 


0.00 


0.90 


0.90 


1.00 


9.09 


16,36 


9.56 


0.47 


0.90 


0 . 85 


10 . 0 


24.37 


17.53 


23.94 


- 0.43 


0.90 


0.85 


6.92 


21.16 


17.18 


21.21 


- 0 .05 


0.90 


0 . 85 


5.00 


18.42 


16.64 


19.06 


0.64 


0.90 


0 .85 


2.00 


13 . 50 


16.49 


14.10 


0.60 


0 . 90 


0.85 


1.00 


11.31 


17.35 


11.22 


- 0.09 
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Table 4 ( continued ) 



l/p 2 


C^/Cj 




®DW 


K 


®DC 




0.80 


0.95 


10.0 


16.33 


17.54 


16.03 


- 0 .30 


0 . 80 


0.95 


6.92 


14.20 


17 . 21 


14.20 


0.00 


0 . 80 


0.95 


5 .00 


12.79 


17.25 


12.76 


- 0.03 


0 . 80 


0 .95 


2 .00 


9.82 


17.91 


9 . 44 


- 0 .38 


0.80 


0.95 


1.00 


7.58 


17.36 


7.52 


- 0.06 


0 . 80 


0 .90 


10 . 0 


20.46 


16.76 


21.02 


0.00 


0.80 


0.90 


6.92 


19.14 


17.70 


18.62 


- 0 . 52 


0.80 


0 . 90 


5.00 


16.33 


16.80 


16 . 54 


0 . 21 


0.80 


0.90 


2 . 00 


12.02 


16.72 


12.38 


0.36 


0.80 


0 .90 


1.00 


9.82 


17.16 


9.86 


0.04 


0.80 


0.85 


10 . 0 


25.03 


17.47 


24.67 


- 0.36 


0 . 80 


0.85 


6.92 


22 . 44 


17.68 


21.85 


- 0 . 59 


0.80 


0.85 


5.00 


19.14 


16.78 


19.64 


0 .50 


0.80 


0.85 


2 .00 


14.20 


16.83 


14.53 


0.33 


0.80 


0.85 


1.00 


11.31 


16.84 


11.56 


0.25 


0.70 


0.95 


10 . 0 


17.07 


17.72 


16.59 


- 0.48 


0.70 


0.95 


6.92 


14.20 


16.64 


14.70 


0 . 50 


0.70 


0.95 


5.00 


13 . 50 


17.60 


13.21 


- 0.29 


0.70 


0.95 


2.00 


9 . 82 


17.31 


9.77 


0.05 


0.70 


0.95 


1.00 


8.36 


18.51 


7.78 


- 0 . 58 


0.70 


0.90 


10 . 0 


21.16 


16.75 


21.75 


0 .59 


0.70 


0.90 


6.92 


19.80 


17.70 


19.27 


- 0 .53 


0 . 70 


0.90 


5.00 


17.07 


16.98 


17.31 


0 . 24 


0.70 


0.90 


2.00 


12.79 


17.20 


12.81 


0.02 


0.70 


0 . 90 


1.00 


10.59 


17.88 


10.20 


- 0.39 


0.70 


0.85 


10.0 


25.64 


17 .30 


25.52 


- 0.12 


0 . 70 


0.85 


6.92 


23.12 


17.61 


22.61 


- 0.51 


0 . 70 


0 . 85 


5.00 


19.80 


16.78 


20.32 


0.52 


0 . 70 


0 . 85 


2 . 00 


14.95 


17.13 


15.03 


0.08 


0 . 70 


0 . 85 


1.00 


12.02 


17.30 


11.96 


0.06 
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Table 4 (continued) 



l/p2 


C 1 /C 2 




®DW 


K 


®DC 


46, 


0.50 


0.95 


10.0 


17.74 


16.91 


18.07 


0 


0.50 


0.95 


6.92 


15 . 64 


16.82 


16.01 


0 


0.50 


0.95 


5.00 


14.95 


17.90 


14.38 


-0 


0.50 


0.95 


2.00 


9.82 


15.89 


10.64 


0 


0.50 


0.95 


1.00 


9.09 


18.48 


8.47 


-0 


0.50 


0.90 


10.0 


22.44 


16.31 


23.69 


1 


0.50 


0.90 


6.92 


21.80 


17.89 


20.99 


-0 


0.50 


0.90 


5.00 


18.42 


16 . 82 


18 .86 


0 


0.50 


0.90 


2.00 


14.20 


17.53 


13.95 


-0 


0.50 


0.90 


1.00 


11.31 


17 . 54 


11.10 


-0 


0.50 


0.85 


10.0 


27.20 


16.85 


27.80 


0 


0.50 


0.85 


6.92 


25.03 


17.50 


24.62 


-0 


0.50 


0.85 


5.00 


21.80 


16.96 


22.13 


0 


0.50 


0.85 


2.00 


17.07 


17 .96 


16.37 


-0 


0.50 


0.85 


1.00 


13.50 


17.84 


13.03 


-0 










Average 


K = 17.22 








Standard (N-1) 


deviation 


= 0.51 










(N-1) 


variance 


= 0.26 





.33 

.37 

.57 

.82 

.62 

.25 

.81 

.44 

.25 

.21 

.60 

.41 

.33 

.70 

.42 
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APPENDIX D 



Case 1 

2.12 



Measured Measured 



10 



100 



max 


a 


®DW 


BW 

w 


9d 


BW 


1 


1.60 


24.5 


10.33 






1 


1.10 


22 . 0 


11.34 






1 


0.83 


19.0 


12.88 






1 


0 . 66 


16.5 


15.27 


13 . 5 


15 


1 


0.55 


14.5 


18.38 


15 . 5 


15 


4 


1.60 


23.0 


6.28 






4 


1.10 


22.1 


8.19 






4 


0.83 


19 . 0 


12.20 






4 


0 . 66 


16.0 


13.88 


13.5 


15 


4 


0.55 


15.0 


14.71 


15.5 


15 


1 


1.60 


17.5 


19.26 






1 


1.10 


16 . 0 


20.29 






1 


0.83 


16.5 


20.90 






1 


0.66 


13.5 


22.41 


13.5 


15 


1 


0.55 


12 . 5 


23.19 


15.5 


15 


4 


1.60 


16.0 


13.85 






4 


1.10 


13.5 


14.96 






4 


0 . 83 


13.0 


15.36 






4 


0.66 


12.0 


15.67 


13.5 


15 


4 


0.55 


12.0 


16.33 


15.5 


15 


1 


1.60 


16.5 


22 . 68 






1 


1.10 


15.5 


23.27 






1 


0.83 


14.0 


23.91 






1 


0.66 


13.0 


24.81 


13 . 5 


15 


1 


0.55 


11.5 


26 . 00 


15.5 


15 


4 


1.60 


12.0 


15.18 






4 


1.10 


12.0 


14.93 






4 


0.83 


11.5 


15.25 






4 


0.66 


11.0 


15.94 


13 . 5 


15 


4 


0.55 


10.5 


16.95 


15.5 


15 


C 1 /C 2 = 


0.8854 


P^/P 


2 = 


8415 


6 = 1.52 
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Case 2 



X 

s 



3.50 



10 



100 







®DW 


BW 

w 


Measured 


Measured 


max 


a 




BW 


1 


1.21 


21.5 


18 . 31 






1 


0.83 


17.5 


20 . 60 






1 


0.63 


14.5 


24.19 






1 


0.50 


12.5 


26.62 


13 . 3 


15.5 


1 


0.42 


10.0 


27.41 


14 . 0 


14.38 


4 


1.21 


22.5 


15.04 






4 


0.83 


16.0 


19.24 






4 


0.63 


14.0 


19.32 






4 


0.50 


13.0 


20.34 


13.3 


15.5 


4 


0.42 


11.5 


22.75 


14.0 


14.38 


1 


1.21 


18.0 


24.19 






1 


0.83 


14.5 


25.13 






1 


0 .63 


14.5 


28.69 






1 


0.50 


12.0 


29.56 


13.3 


15.5 


1 


0.42 


9.5 


30.18 


14.0 


14.38 


4 


1.21 


16.0 


16.07 






4 


0.83 


14.0 


17.97 






4 


0.63 


14.0 


20.43 






4 


0.50 


12 . 5 


22.70 


13.3 


15.5 


4 


0.42 


11.5 


24.96 


14.0 


14.38 


1 


1.21 


18.5 


26 . 64 






1 


0.83 


16.0 


27 . 72 






1 


0.63 


14.0 


29.29 




15.5 


1 


0.50 


11 . 5 


30 . 89 


13.3 


1 


0.42 


9.5 


31.55 


14.0 


14.38 


4 


1.21 


14.0 


17 . 65 






4 


0.83 


13.5 


17.94 






4 


0.63 


12.5 


19.37 




15.5 


4 


0 . 50 


11.5 


21.69 


13.3 


4 


0.42 


11.0 


23.80 


14.0 


14.38 


Co = 


0.8854 


Pl/P, 


= 0. 

1 


8415 6 = 2.51“’ 
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Case 3 



6.15 



10 



100 





a 


®DW 


BW 

w 


Measured 

®D 


Measured 

BW 


1 


0,58 


13 . 5 


30.13 






1 


0,40 


8.0 


32.50 






1 


0. 30 


7.0 


33.47 






1 


0.24 


7.0 


35.75 


16.4 


25.5 


1 


0,20 


7.5 


38.56 


20 


21 


4 


0,58 


14.0 


24.62 






4 


0.40 


10.5 


30.65 






4 


0.30 


10,5 


35.06 






4 


0.24 


10.5 


41.51 


16.4 


25.5 


4 


0.20 


11.0 




20 


21 


1 


0.58 


14,5 


24.60 






1 


0.40 


12,5 


30.47 






1 


0.30 


11.5 


35.89 






1 


0.24 


11.5 


41.62 


16 . 4 


25.5 


1 


0,20 


12.0 




20 


21 


4 


0,58 


14.5 


33.94 






4 


0.40 


9.0 


35.10 






4 


0,30 


7.5 


36.04 






4 


0.24 


7.5 


38 . 04 


16.4 


25.5 


4 


0.20 


7.5 


40.78 


21 


21 


1 


0.58 


14.0 


35.19 






1 


0.40 


8.5 


35.83 






1 


0.30 


7.0 


36.47 






1 


0.24 


7.0 


38 . 37 


16.4 


25.5 


1 


0.20 


7.5 


40.82 


20 


21 


4 


0,58 


13.5 


25.14 






4 


0.40 


11.0 


30.32 






4 


0.30 


10.0 


35.38 






4 


0,24 


10.0 


40.97 


16.4 


25.5 


4 


0.20 


11.0 


35.02 


20 


21 


^/C2 = 0, 


8854 


P^/P 


2 = 0- 


8415 B 


= 4.41“’ 
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Case 4 



%ax 


a 


DW 


1 


1.31 


18. 


1 


0.90 


14. 


1 


0.68 


13 . 


1 


0.54 


10. 


1 


0.45 


8. 


4 


1.31 


17 . 


4 


0.90 


15. 


4 


0.68 


14. 


4 


0.54 


12. 


4 


0.45 


11. 


1 


1.-21 


17. 


1 


0.90 


14. 


1 


0.68 


12. 


1 


0.54 


10. 


1 


0.45 


8. 


4 


1.31 


14. 


4 


0.90 


13. 


4 


0.68 


12 . 


4 


0.54 


11. 


4 


0.45 


10. 


1 


1.31 


16. 


1 


0.90 


13 . 


1 


0.68 


11. 


1 


0.54 


9. 


1 


0.45 


7. 


4 


1.31 


12. 


4 


0.90 


11. 


4 


0.68 


11. 


4 


0.54 


10. 


4 


0.45 


9. 



C^/C2 = 0.9214 



BW 

w 


Measured 

9d 


Measured 

BW 


15.88 

17.62 

20.75 

22.84 


12.67 


12.88 


24.58 


14.7 


13.5 


9 .65 
13 .38 
17.32 
21.72 


12.67 


12.88 


25.21 


14.7 


13.5 


21.32 

22.55 

24.35 

26.14 


12.67 


12.88 


26.62 


14.7 


13.5 


13.80 

15.71 

17.96 

20.56 


12.67 


12.88 


22.82 


14.7 


13.5 


22.58 

23.71 

25.18 

26.56 


12.67 


12.88 


26.76 


14.7 


13.5 


15.05 

15.87 

17.27 

19.32 


12.67 


12.88 


21.15 


14.7 


13.5 


= 0.' 


9073 


6 = 2.66° 



5 

5 

0 

0 

0 

5 

5 

0 

0 

0 

0 

0 

0 

0 

0 

5 

0 

0 

5 

0 

0 

5 

5 

5 

5 

5 

5 

0 

0 

0 

/P2 
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FUNCTION OF ANGLE CAPITAL THETA (THETAK) 
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FROM THE BOTTOM DOMAIN CALCULATIONS 
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